Crumpling of a stiff tethered membrane 
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A first-principles numerical simulation model for crumpling of a stiff tethered membrane is introduced. In our 
model membranes, wrinkles, ridge formation, ridge collapse, as well as the initiation of stiffness divergence, are 
observed. The ratio of the amplitude and wave length of the wrinkles, and the scaling exponent of the stiffness 
divergence, are consistent with both theory and experiment. We observe that close to the stiffness divergence 
there appears a crossover beyond which the elastic behavior of a tethered membrane becomes similar to that 
of dry granular media. This suggests that ridge formation in membranes and force-chain network formation in 
granular packings are different manifestations of a single phenomenon. 
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Crumpling of, e.g., sheets of paper, is an everyday phe- 
nomenon that one easily passes without much extra attention. 
For physicists, however, crumpling has for a long time been a 
topic of intense research. In addition to fundamental statisti- 
cal mechanic properties QJ, issues such as scaling of strength 
and energy 12 0L geometry and singularities an d even 

acoustic emission |6] and crumpling related phenomena in 
soft matter systems |7|], have become topics of increasing in- 
terest. In this article, we present a first-principles numerical 
simulation model for crumpling of a stiff tethered membrane, 
and relate the problem to formation of force chain networks 
in granular media 

Crumpling of a thin elastic sheet or a stiff membrane de- 
mands deformation energy. At very small strains a membrane 
is uniformly compressed. Thin structures are, however, sen- 
sitive to buckling. At first buckling appears as wrinkles. The 
wave length and the amplitude of the wrinkles depend on the 
conditions of the external loading and on the dimensions of 
the membrane 1 As crumpling proceeds the deforma- 
tion energy begins to concentrate in narrow ridges and con- 
ical peaks Q 0]. Upon crumpling of a piece of paper, for 
example, the strain at the ridges and peaks become so large 
that irreversible plastic yielding takes place. After the paper is 
stretched out again, the ridge pattern can clearly be observed. 

Another characteristic feature of membrane crumpling is 
that the effective stiffness of a crumpled membrane increases 
fast when the degree of crumpling becomes large. The ulti- 
mate limit of crumpling would be to press a membrane into a 
more or less cubic form with a volume corresponding to that 
of the volume of the membrane itself. In practice this limit 
can hardly be reached. Squeezing a piece of paper as hard 
as possible between the hands, results in a ball that typically 
contains about 75% air II ill . The behavior of the effective 
stiffness has been investigated by carefully monitored exper- 
iments, and for a crumpled membrane it displays a power- 
law divergence as the dense packing limit is approached 1 1 1 ]. 
The suggested explanation to the power-law divergence is that 
the total ridge length diverges in a scale invariant way as the 



porosity inside the crumpled membrane vanishes. 

In order to study wrinkles, ridge formation, and the effec- 
tive stiffness in crumpled membranes, we constructed a nu- 
merical model for stiff tethered membranes. The tethered 
membrane consists of frictionless spheres with stiffness Y s 
arranged in a triangular lattice. Each sphere has a mass m 
and moment of inertia /, and is connected to its neighbors by 
massless elastic beams which have both bending and tensile 
stiffness. It is reasonable to use a higher Young's modulus for 
the spheres (Y s ) than for the beams (Y&). In the simulations 
here we chose Y s = 10 7 and Y& = 5 x 10 3 . The membrane 
is placed horizontally (i.e., in the xy plane) in the middle of a 
rigid cube of size X 3 (t), where X is time dependent. Small 
random fluctuations of the z component of the centers of mass 
of the spheres are introduced to avoid perfect symmetry. New- 
tonian dynamics is applied. 

The environment inside the box is set to be strongly dissi- 
pative (the membrane can be thought of being immersed in 
a highly viscous fluid). This means that a force proportional 
to velocity, (dx), slows down the motion of each sphere. The 
membrane dimensions are considered to be large enough for 
temperature fluctuations to be negligible (e.g. a paper sheet). 
For computational reasons the membrane cannot contain more 
than about 10 5 degrees of freedom, which sets a limit on its 
size. Wrinkles and the onset of ridge formation is neverthe- 
less observed. Also the beginning of a power-law divergence 
of the effective stiffness with an exponent close to the one ob- 
served experimentally is seen. However, as porosity reaches 
a low enough value, the effective stiffness begins to deviate 
from the experimental power law. At low porosities the ef- 
fective stiffness of a tethered membranes behaves in the same 
way as that of a packing of solid spheres 1 12]. 

Figure \l\ shows a snapshot from a simulation displaying a 
few wrinkles at t = 200. By minimizing the elastic energy it 
can be demonstrated that the relation between the wave-length 
(A) and the amplitude (A) of the wrinkles is generally given 
by A/X — (A/W) 1 / 2 , where W is the width and A the com- 
pression of the membrane perpendicular to the wrinkles 
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FIG. 1: Snapshots of a crumpling membrane. Top views (above) and 
side views (below) at times t = 200, 500, 700 are displayed. At 
t = 200 wrinkles can be observed. At t = 500 ridges have formed 
and begun to collapse. At t = 700 the membrane is at the beginning 
of the stiffness scaling regime, (X — X c )/X c « 1.6. 




FIG. 2: The ratio of wrinkle amplitude and wave-length (A/ A) as 
a function of membrane compression A. The simulation result is 
compared to the theoretical prediction Aj A oc vA. 



FIG. 3: Contour maps (at t = 400) of the spatial locations of (A) the 
highest local bendings p(x, t) < 0.6, (B) the highest elastic energy 
of the beams E^(x, t), and (C) the highest contact elastic energies of 
the spheres E s (x, t) . 
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FIG. 4: (A) The total force F on the box surrounding a membrane 
scaled by sphere Young's modulus Y 3 and linear membrane length 
W as function of (X — X c )/X c . Data is shown for membranes of 
linear size W = 40, 80. The solid line shows the experimentally 
observed Q power-law divergence F(X) ~ (X - A" c ) -1 ' 85 . (B) 
F/(WY S ) as function of solid volume fraction <j>. The scaling of 
membrane stiffness (cj> c — 0) -1 ' 85 , with (j> c w 0.76, and the scaling 
of granular stiffness (<f> — </> c ) 1,62 , with cf> c ~ 0.67 are compared to 
the simulation results. 



This relation is tested in Fig.[2]for the model membranes. The 
square-root relation between Aj A and A holds up to A w 6 
(here W « 23.4). 

Deformation of ridges in controlled geometries has been 
thor ough ly inve stiga ted both numerically and theoretically |2, 
0, S E3- E3- EH E3l • The elastic deformation energy of a 
ridge of length L scales as ^(L/S) 1 / 3 , where 5 is the thick- 
ness of the sheet and k its bending modulus. By taking advan- 
tage of this result, and following the derivation in Ref. full , 
the scaling of stiffness of a membrane can be estimated: A 
crumpled sheet is divided into facets by the ridges. Ridges 
surrounding a facet of size I? have a length proportional to L. 
Such facets should fill a space of size L 3 . This leads to a bulk 
stiffness divergence of the form K — —VdP/dV oc V^ 11 ^ 3 , 
when the volume V surrounding the sheet decreases under the 
influence of pressure P. This is not quite consistent with ex- 
perimental findings. In Ref. full K was reported to scale as 
K ex (V c - V)~ a with a sa 2.85. Notice that this scaling rela- 
tion can also be expressed in the forms F oc {4> c — <fi)~ a+1 and 
F oc (X — X c )~ a+1 , where F is the external force applied 
on the membrane, <fi the solid volume fraction, and X(t) the 
time dependent linear dimension of the crumpling membrane. 

To investigate possible ridge formation in the simulations, 
we compare the local bending to the local energy. We define 
the sum of the dot products of nearest neighbor bonds along 



the three principal directions in a triangular lattice as a mea- 
sure of the local bending of the membrane, 

6 i=l L i l i 

This is then compared to the local elastic energies of the 
beams 

E b {x,t)= V / F{Y b ,x(t))dx(t). (2) 

In Eq. Q, l m and l n are the lengths of the beams on the 
opposite sides of a lattice site. In Eq. (0, the elastic force 
F(Y"b, x(t)) is integrated over the corresponding displacement 
and summed over the six degrees of freedom for each site. 

In addition, the elastic deformation energy of spheres in 
contact can be used to identify ridges. It can be written as 

n 1 

E s (x,t) = J2^s, (3) 

i=l 

where Y s is the stiffness of the spheres and 5 S the virtual 'over- 
lap' of spheres in contact. The sum is taken over all contacts. 

Figure|3]shows contour plots of the lowest values of p(x, t) 
and the highest values of Ef,(x,i) and E s (x,t) at an early 
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stage of the crumpling process (t = 400). It is quite obvi- 
ous that Eb and E s reach their highest values at the locations 
of high local bending. The patterns are quite complex, but 
elongated structures, i.e. ridges, can be seen in the figure. As 
crumpling proceeds, ridges grow and collapse and the mem- 
brane folds. For later crumpling stages it is no longer possible 
to visibly detect ridges as the patterns become too complex. 

The effective bulk stiffness of a membrane can be moni- 
tored, e.g., through its stress-strain curve. This is simply the 
total force F on the walls of the box surrounding the mem- 
brane, divided by the linear size W of the membrane (F/W) 
as a function of X(t). Initially F(t)/W increases fairly lin- 
early until the membrane begins to buckle and fold. As the 
dense packing limit is approached, F/W begins to increase 
rapidly. Finally, when the membrane almost fills the box, a 
saturation regime is entered. 

It is important to notice that there are three types of dense 
packing limits: 1) The tightest packing is achieved when the 
membrane volume is equal to the volume of its container. In 
this limit, the solid volume fraction is <f>i = 1.0. 2) Hard- 
core spheres. Here, the tightest packing is an FCC lattice with 
<t>2 = 7i"v2/6 w 0.74. 3) A random dense packing of spheres 
with 03 w 0.63. For membranes it seems that the stiffness 
divergence takes place at <f> c w 0.75 which is close to <j>2- The 
regime of saturating stiffness is, however, entered well before 
this limit is reached. 

Figure 0\ shows F/(Y S W) for membranes with 10 4 and 
4 x 10 4 degrees of freedom (each simulation demands about 
10 6 — 10 7 time steps). The simulation results are compared 
to the experimental result reported by Matan et al. lUlL who 
found that F{X) oc (X - X c )- a+1 , with a « 2.85. Fig- 
ure |4j\ shows good agreement between simulations and ex- 
periment until saturation sets in at (X — X c )/X c w 0.03. 

At saturation there is a crossover to a different type of stiff- 
ness behavior. Beyond saturation the beams no longer affect 
the stiffness of the membranes. Instead, the membranes begin 
to behave like granular packings lll2ll . The stiffness of a 3D 
packing of solid spherical grains scales like (<fi — 4> c ) , with 
/3 « 1.62, and the critical volume is <\> c « 0.6 — 0.7 depend- 
ing on the form of interaction between the grains. Figure @J3 
shows a comparison between simulations and the scaling be- 
havior of F(4>) in the two mentioned regimes. 

By following the energy contents in a membrane from 
the initially flat membrane to a fully crumpled state, it is 
possible to determine the energy component responsible for 
the stiffness divergence seen in Fig. |4] The external energy 
J t F(t) \X(t) — X({))] dt is compared to the internal energy 
of the membrane in Fig.|5J\. There is no difference between 
the external and internal energies as energy conservation also 
demands. 

The internal energy can be split into five terms that are re- 
lated to the force terms in the general equation of motion 

mx + dx + Kx + Y S 5 S = f ext . (4) 

Here, x is the displacement, m the mass, d the damping co- 
efficient, K the stiffness matrix, S s the deformations (or vir- 
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FIG. 5: (A) Internal and external energies as functions of time t. (B) 
E 3 and E ext as functions of t. (C) Ed, Et and E m as functions of t. 

tual overlap) of spheres in contact, and f ext are forces applied 
on the walls of the box surrounding the membrane. The re- 
lated energy terms are the following: (1) The elastic energy 
of the beams, Eb. (2) The compression energy at the con- 
tacts between spheres, E s . (3) The compression energy of the 
spheres against the external walls, E ex t — 0.5^2Y W S^,. (4) 
The energy dissipated through damping, Ed = J t Fddx. (5) 

The kinetic energy of the spheres, E rn = 0.5 m ^ ■ These 
energy components are shown separately as functions of time 
in Figs.[5j} and|5jT. The diverging terms are E s and E ext . The 
divergence of E ext is rather trivial as it is just the deformation 
energy of the box surrounding the membrane. It is thus the 
energy E s that causes the stiffness divergence. This explains 
to some extent why there is a crossover to the stiffness behav- 
ior of granular packings. When the deformation energy at the 
contacts between spheres dominate over other energy terms, 
the beams no longer play a role in the membrane stiffness. 

Since there are two different scaling behaviors of stiffness, 
and both regimes are dominated by E s , one would expect that 
at the crossover point there appears a qualitative change in the 
function describing the distribution of contact energies of the 
spheres. There is, however, no indication of such a qualita- 
tive change as is demonstrated in Fig. [6] which means that, at 
the crossover, ridges are simultaneously force chains. When 
comparing the topology of the network of ridges in crumpled 
membranes and that of the force-chain network in dry granu- 
lar materials |3|9[, one notices that both are essentially one- 
dimensional lines of localized deformation energy, and form 
complicated networks in three-dimensional space. Based on 
these observations, it is plausible to draw the conclusion that 
the ridge formation and force-chain formation are not sepa- 
rate phenomena, but different manifestations of a single phe- 
nomenon. The crossover between the two then only marks 
the region where the membrane thickness begins to limit the 
density of network lines. 

In conclusion, we have demonstrated for the first time that 



4 




Es/<E^> 4> 

FIG. 6: (A) Distribution functions for the deformation energies at 
contacts between spheres for different solid volume fractions <f> — 
0.04, 0.37, 0.57, 0.76. The distribution functions are compared to 
the equilibrium distribution exp(— E s ) 1 18]. (A) E s as a function of 
cj>. The crossover that takes place at (j> ~ 0.72 is not visible in the 
data. 

a numerical model of a tethered membrane can reproduce the 
theoretically predicted wrinkling and ridge formation, and ex- 
perimentally observed stiffness divergence. The early-stage 
ridges can be traced both by the elastic deformations of the 
beams and by the contact deformations of spheres in a teth- 
ered membrane. The divergence of the effective stiffness is a 
result of the divergence of the contact energy of the spheres. 
This eventually leads to a crossover in the stiffness between 
two different scaling behaviors. When the solid volume frac- 
tion is far below that of the random dense limit of granular 
packings, the scaling of stiffness is 'membranic', and when 
the random dense packing is approached, the stiffness scal- 
ing becomes 'granular'. There is no obvious difference in the 
distribution of deformation energy in the two regimes. It thus 
seems that deformation ridges in membranes and force chains 
in granular packings are different manifestations of the same 
phenomenon. 
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